Casimir effect in the nonequilibrium steady-state of a quantum spin chain 
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We present a fully microscopies-based calculation of the Casimir effect in a nonequilibrium system, 
namely an energy flux driven quantum XX chain. The force between the walls (transverse-field 
impurities) is calculated in a nonequilibrium steady state which is prepared by letting the system 
evolve from an initial state with the two halves of the chain prepared at equilibrium at different 
temperatures. The steady state emerging in the large-time limit is homogeneous but carries an 
energy flux. The Casimir force in this nonequilibrium state is calculated analytically in the limit 
when the transverse fields are small. We find that the the Casimir force range is reduced compared 
to the equilibrium case, and suggest that the reason for this is the reduction of fluctuations in the 
flux carrying steady state. 
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I. INTRODUCTION 



a. 

Historically, the Casimir force [l[ is the effective interaction that develops between two ideal conductors in the 
vacuum due to the quantum fluctuations of the electromagetic field. The zero-temperature case was soon generalized 
by Lifschitz Q who also calculated the effective force induced by thermal fluctuations. Since then, these theoretical 
predictions have received, as decades elapsed, qualitative and then quantitative confirmations [H, [ij. For the electro- 
magnetic field, the thermal component of the Casimir force is orders of magnitude weaker than its zero-temperature 
counterpart. This accounts for the relatively late confirmation of the thermal Casimir force [H, H|. The Casimir 
force has also appeared in other frameworks, such as low dimensional quantum liquids It was found there, for 

' example, that the Casimir force between magnetic impurities displays Friedel-like oscillations. This large body of the 
, existing literature is mainly devoted to equilibrium situations. 

In this paper, we wish to investigate the interplay of the long range effective interactions produced by an energy flux 
■ * \ running through the chain and the fluctuation mediated interactions between two point-like defects located a distance 
I apart. Following the accepted terminology, the effective interaction between the defects will be called as Casimir 
force. We have been inspired by a recent series of works by Antezza et al. [If! [l4| in which a similar problem was 
t-H . considered, namely the calculation of the Casimir force produced by the electromagnetic field in between two parallel 
plates kept at unequal temperatures. This work followed earlier investigations by Polder and Van Hove [ll|, recently 
. i \ reviewed by Volokitin and Persson [l2| and Dorofeyev [Hj]. However, in contrast to [13, [HI, m our study there will 
be no phenomenological input in the form of linear response coefficients; it will be fully microscopicallly based, both 
for modeling the walls and the bulk of the system. Our theoretical laboratory will be an integrable spin chain(the 
XX quantum chain) where similar calculations have been carried out before for purely equilibrium situations [a [7 , |9 1 . 
The results we have obtained notably differ from those obtained in a different framework by Antezza et al. UA | ■ 
In particular, we find that in constrast to their work, even in the simplest case of two parallel plates, in the presence 
of an energy flux, the Casimir force does not appear to be the average of the contributions of two equilibrium Casimir 
forces, corresponding to the two imposed temperatures. We, on the other hand, do not have a simple decomposition 
due to the fact that the energy flux is the quantity which governs the value of the force, and in our case it clearly 
emerges that the presence of the flux decreases the magnitude of the force. 

Our results are presented in the following order. Section |H] briefly reviews the XX chain and describes the properties 
of nonequilibrium steady-states arising from preparing the system with a step-like temperature profile. In Section [III1 
we introduce two magnetic impurities and we show how the steady-state properties are modified in their presence. 
In Section IIV1 we derive the Casimir force exerted by one impurity upon the other and comparison with earlier 
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phcnomenological calculations are presented. Conclusions and future research directions follow in Section [Vj Finally, 
in appendix lAl we present a physicist's derivation of the nonequilibrium steady-state properties, previously obtained 
via more complex C* algebraic methods. 



II. TRANSVERSE XX CHAIN 



A. Equilibrium 

The transverse XX chain is one of the simplest quantum systems displaying long-range correlations and the asso- 
ciated large fluctuations. It is defined by the Hamiltonian 

i 

where Si — \di and af (a — x,y,z) denote the three Pauli matrices at sites — N, — 1, 0, N — > oo of a d = 1 



chain, and h is the transverse field in units of the coupling J and J = 1 is set in the following. The standard way 15 1 
to study the transverse XX chain is to resort to the Jordan- Wigner transformation which maps the spin chain onto a 
one-dimensional system of free fermions with energy spectrum e q = — cos q — h: 

H = Y,^A c <i ( 2 ) 

<? 

where c q is the Fourier transform of c- L = o~ ^rij<i-i( — with the wavenumbers in the range — it < q < it. 
As far as large-distance or transport properties are concerned, most of the relevant large scale physics is governed 
by the modes q in the vicinity of the Fermi level ±k determined from cosk = —h. We shall thus often resort to 
the approximation of effective relativistic fermions which consists in linearizing the dispersion relation around ±k (by 
setting q = ±n + k) and considering the modes with q > or q < as two independent families of relativistic fermions, 
namely the left and the right movers, with velocity c = sin k. Further details on the validity of such a description can 
be found e.g. in [l6[. A phenomenological cut-off A is imposed, when needed, on the new k modes. 



B. Steady state with energy flux 

Our goal is to investigate the nonequilibrium states of H which carry a given energy flux (Je) 0. Here the energy 
flux Je is given by 

Je = / j sin qeqCgCq (3) 

One way to achieve a current carrying steady state has recently been discussed by Ogata [1 71 ] , who built upon earlier 
studies by Araki [l8| , Tasaki , and Pillet and Aschbacher [ijj] . Initially, the chain is prepared in the following way: 
its left side, say from — oo up to site j = is in thermal equilibrium at inverse temperature while its right hand 
side from j = 1 up to +00 is in equilibrium at inverse temperature /?2- Both halves are initially disconnected. Then 
at t = contact is made through the (j = 0,j = 1) bond and the system eventually settles into a nonequilibrium 
steady-state. Due to its infinite thermal conductivity, the temperature profile is flat in the central region whose extent 
expands at finite velocity. We shall not be concerned here with the dynamics of the formation of the central region and 
with the front propagation issues (2]H23 |. We focus on the asymptotic homogeneous steady-state of the expanding 
central region of the system, where the fermion occupation number in the steady-state is given by 

/ t v 6{k) 0(-k) 

The above result can be derived mathematically rigorously [17] . a physicist's proof of (@| is provided in appendix[XJ It 
will often prove convenient to introduce (3 = £i±* and e' k = 6(k)^-ek + 0(—k)^-ek, so that Fk appears as the effective 
Fermi-Dirac occupation number at temperature f3 of free fermions with energy spectrum e' fe . The discontinuity of e' k 
is the consequence of the long-range effective interactions in the steady-state, as discussed in [T3|. Ogata has also 
discussed the dependence of (Je) on h, fii and fi 2 in light of the experimental literature (27|. We note that there exist 
other ways to produce a nonequilibrium steady-state (see [U [2^, [2(| ) , but we shall not consider these here [28j ] . 
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III. A SPIN CHAIN WITH TWO IMPURITIES 



A. About the force 



In our case, the walls of the Casimir setup will be two magnetic impurities at lattice sites z f£/2 with strength 5h\/ 2 
which enter the Hamiltonian through additional terms of the form 

SH = -8h 1 c i _ e/2 c_ e /2 - Sh 2 cl /2 c e / 2 (5) 

The task is to determine the effective force between those impurities. Since we are ultimately interested in a nonequilib- 
rium setting, we must circumvent the methods used beforehand in 0, 0, [H| that relied on determining a ground-state 
energy or a free energy. We adopt the following definition of the force on a lattice. Let Hj be the energy density at 
site j, 

Hj = -~{c]c j+1 + c] +1 cj) - hc]c,j (6) 
The force F felt by the defect located at site — £/2 is given by 

F = ~{H_ i/2 -i - H_ e/2+ i) (7) 

where the average is over the nonequilibrium steady-state. The local energy being quadratic in the fermionic operators, 
it is clear that in order to calculate F, it is sufficient to know the two-point equal time Green's function Qij{t,t') = 
(ci(t)Cj(t)} as calculated in the next Subsection. Note that there is some arbitrariness in defining a force directly on 
the lattice. In a continuum limit, e.g. when one focuses on wave vectors close to the Fermi level, such differences 
become irrelevant, and our definition matches the expression of the force one would obtain directly from a continuum 
theory. 



B. Green's functions 

In order to determine the Green's function Qij(t,t') = (cj(t)ct(t')), we shall need another Green's function Gij(t,t') 
introduced through the Heisenberg picture, 

c,(t) = ]TG m (t,0)c„(0), c}(t) = Y,G* mj (t, 0)4(0) (8) 

n m 

which leads to 

SaM = G in (t,0)G* mj (t', 0)^(0)4(0)) (9) 

In the latter formula, we take as the initial state the stationary-state. The Green's function G has to be calculated in 
the presence of defects, which can be carried out by a variety of methods, one of them being presented in appendix FBI 
The results can most conveniently be expressed in terms of the Green's function <7 mn (oj) in the absence of defects, 

5m«H = / ~z : : (10) 

J 2ir —iuj + e q 

where uj is conjugate to it. We find that the Fourier transform of Gij reads, in the weak magnetic defect limit, 

Gij =gij + Shigi^g-j + 8h 2 gi+g+j 

+ dhfg^g^jg — + 5h 2 2 g l+ g +j g ++ (11) 
+ 5hi8h 2 (g l+ g^ :j g + - + g^g +j g- + ) + 0(Sh 3 ) 

where ± is for ±4. The alternative limit Shi, Sh 2 oo would be corresponding to freezing the degrees of freedom of 
the defects, which would bring us closer to the original setting of Casimir. In terms of space Fourier transforms, this 
expression becomes 



9q9q' 



r r\k ( 12 ) 
+ / — (W^'-^ +Sh 2 e-^'-V^ (dh^-^i +5h 2 e~ 1 ^^ g q g q ,g k + 0(Sh 3 ) 
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We now use that the steady-state Green's function Q is given in (j9)), where the initial state (c n (0)cj n (0)) is the 
current-carrying steady-state obtained from a step temperature profile at t = — oo and G includes the presence of 
weak impurities as given in In the infinite time limit, we show in appendix IB1 that 



G(q,q',t) - F Q$ q ,q' - 

dk 



(13) 



1 



Fn - F k 



F k 



E q — Ek £q> — £k 

Since we are interested only in the interactions between the defects, only the cross Shi x 8h 2 term contributes to 
Gx(q,q') when calculating the corresponding force: 



G x (q,q') =26h 1 5hi 



dfc cos(q + q' - 2k) { 
2ir e„ - e„> 



It is useful to introduce the function 7 (r, u) = J ^e lqr j q (Lo) whose Fourier transform is defined by 



F q -F k 



F q ' — F k 



£q — £k 



£«' — £k 



With this notation, we have 

G x {q,q') =25hx8h 2 
In real space, this becomes 



7 9 (w) 



du> dk 
2tt 2tt ' 



— vjJ + e' 



(q + q' -2k)- 



Jqi^hq' M 7 fc(w) 



(14) 



(15) 



(16) 



Gx(x,y) =6hi5hz 



did 

2^ 



h(£/2 - x, w) 7 (*/2 + y, u)y(-i, u) + j{-£/2 - x, w) 7 (-*/2 + y, w) 7 (*, w)] (17) 



and this form will be used in the calculation of the force below. 



IV. CASIMIR EFFECT: THE EFFECTIVE FORCE BETWEEN THE IMPURITIES 



A. Formal expression of the Force 

In this section we derive the expression of the force in the limit where the impurities are a large distance £ apart. 
In this limit, many of the expressions encountered above simplify significantly. Most notably, the function 7 (r, ui), 
which enters the final expression of the Green's function, reads 



7 (r,w) = e ,;Kr 7l (r,w) +e 4Kr 72 (r,w) 



where we used the following definitions 



and 



„-|ru)|/ci 

ji{r,u) = i (0(r)0(w) - 0(-r)0(-w)) 

ci 



e -\rui\/c2 

72 (r, u) = % (-0(r)0(-w) + 9(-r)B(u)) 

C2 



(18) 
(19) 

(20) 



with c\ — ^-sinK and c 2 = ^sinz-c. In equation (|18l) . the right (1) and left (2) moving fermions account for the C\ 
and c 2 dependent terms, respectively. In terms of 7l (r, w) and 72 (r, w), the force can be written as 

t + 

-<- — y — 

I C\ 



F = 



2/3 



ic 2 
2f3 



J2(-d y + d^e-^-^T^x^,^) 



x=y= 



x=y= 



(21) 
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where the frequency sum is over the ui = 2 "^ l " 1 7r, and where 

ri(a;, y, u) = 5ht6h 2 {e~ 2M l2 {l, wfrxfo - 1/2, - */2, u) + e 2M j 2 (-£, u^y + 1/2, w) 7l (-i + £/2, w)) 

(22) 

and 

r 2 (z, y, w) = Shi 5h 2 (e 24Ki 7i(^, - t/2, uj) 12 (-x - 1/2, u) + e^SlM, u) l2 {y + 1/2, u) l2 (-x + 1/2, u)) 

(23) 

Finally, in the limit f > 1, the force between impurities given by Eq. ([7]) is the discontinuity of the energy density 
across the defect, and in terms of 7i(r, w) and j 2 (r,oj), it can be written as 

F = 2^Ti(*,w)7 2 (-^a;) 

w (24) 
+ >> - + - 7i 72 M,w 

Note that the force can be written in terms of the unperturbed Green's functions of the left and right moving fermions. 
The expression of the force (|2"4"|) is the central result of this work. We now specify this result, first to a known situation 
to make contact with existing results, and second to the physically more interesting case of a current-carrying chain. 



B. Equilibrium (recovering existing results) 

In Q the authors studied the Casimir forces between defects in equilibrium one-dimensional quantum liquids at 
zero temperature, in order to reproduce their results we have to take the (3\ — /3 2 — ft with j3 — > oo in equation (|24l) . 
In this limit the sum over uj can be changed by an integral 4 J2uj ~* I I^j an< ^ a ft er some algebra we find 

F = - iklih J^^ + ^L) (25) 

\ tti sin k 2irz* sin k ) 

which is indeed the interaction force associated to the interaction potential between impurities given by equation (18) 
of reference Q ■ We conclude that in the case of zero temperature the leading term of the interaction force decays as 
l/l and it oscillates with wavelength tt/k. 



C. Out-of-equilibrium, with a heat flux 



In the case /3i ^ /3 2 a heat flux drives the spin chain into a nonequilibrium steady state. The expression for the 
force reads 

45hi5h 2 K sin(2n\e\) ( 1 \ 

F — — 



1ClC 2 



e e - e & 



4irShi5h 2 pcos(2Kz) e 

P 2 C\C 2 



[e * - e f> 

where p — — ^ + ^ and we recall that c\ = 4p sinK and c 2 — fy sin ft. As one can see force decays exponentially with 
a caracteristic length £ = This is in contrast with the zero temperature case where the decay is algebraic. The 
oscillatory factors sin(2«; |£|) and cos(2k€) are the same as in the zero temperature limit. As expected, the /3\ = fi 2 = /3 
thermal equilibrium limit is consistent with previously found results [7] . 



V. CONCLUSIONS 



We have fully determined the Casimir force between two magnetic field defects of weak amplitude in the nonequi- 
librium steady-state of the XX spin chain carrying an energy flux. The overall qualitative behavior is similar to that 
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obtained from equilibrium thermal fluctuations: the Casimir force decays exponentially with the distance in between 
the impurities. However, that decay is notably sharper in the presence of an energy flux than without. Thus we 
conclude that the presence of the energy flux tends to weaken the Casimir force. This can be seen from the cor- 
relation length £ = since, as the system approaches equilibrium fa — > fa, Coutofcq — £cq = — c ^g~^ 2 ' 1 , and thus 
£out of cq < £cq- This strengthens the general picture (30l. l3lj that nonequilibrium fluctuations lead to stiffer systems, 
which, as the present calculation reveals, do not favor fluctuation mediated interactions. 

It is expected that the leading behavior is different in the presence of strong impurities. And so it would be 
interesting to push our investigations further to probe the similarities and the differences with the calculation of 
Antezza et al. [10] bearing on the electromagnetic field between two plates thermalized at unequal temperatures. The 
conceptual issue here is whether different "ensembles" with either a fixed temperature difference or a fixed energy 
flux should lead to the same physical results in a nonequilibrium setting. Another issue of interest is related to our 
spin chain having an infinite conductivity. The question of what happens in realistic systems with nonintegrable 
interactions is an open one which we would like to address in the future. 
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Appendix A: Heat baths at different temperatures 

In her work [17], Ogata resorts to algebraic methods (C* -algebras), to derive the steady-state properties of the 
spin chain whose two semi- infinite halves are initially prepared in thermalized states at T\ and T2 . The goal of the 
appendix is to rederive her results sticking to the standard methods familiar to physicists. 

We start from to semi- infinite XX spin chains defined for n < and n > 1, which are respectively thermalized at 
inverse temperatures fa and fa. The Hamiltonian for the left hand side reads 

= (~^ c « c »+i - \ c l+i c n - hclc n J = ^ dqe q c\c q (Al) 



where the Fourier transform is given by c q = ~iy % J2 n <o sm [(n — l)g]c rl (0 < q < tt) and the energy spectrum is 
e q = — cos q — h. As for the right hand side we also have 

fj{2) = f-lcl iCn+1 - ^4 +1 c„ - hclcA = ^ / dqe q c\c q (A2) 

n>0 ^ /JO 



where the Fourier transform is now given by c q = — i y ~ X)n>i sm [ n 9] c n (0 < Q < tt) and the energy spectrum is of 
course the same. In the initial state of the left hand side of the chain, we have that 

( C t c m ) = - I dqsm[q(n - 1)] sm[q(m - l)]f { q +) , n,m<0 (A3) 
tt Jo 

where fq 1 ^ = c ^ ie 1 g+1 is the Fermi-Dirac occupation number. Similarly, for the right hand side, we have 

(c^c m ) = - I dqsin[qn] sinfgm]/^, n,m > 1 (A4) 
i" Jo 

with fq 2 ^ — c n 2 l q+1 ■ Ogata's result [TtJ states that in the steady-state reached in the large time limit, 

(c{c m ) = ^ ^. e ^n-n)p k (A5) 
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where the occupation number is given by 

1 I 



F k = 



0(k) + ^n^- fc ) ( A6 ) 



In order to prove the result (|A6|) . we start by determining the time dependent Green's function for arbitrary lattice 
sites. Since Ck{t) — e t6kt Ck(0) (the Fourier modes — n < k < it refer to the whole translationally invariant chain), we 
have that 

Cn (t) = P ^e-* kn c k (t) = P ^e^'Te^c^O) (A7) 

J-7T 27T J-n 27r j 

with a similar relationship for cj^t), hence the Green's function evaluated at equal -but finite- times reads 

(4(*)c„(t)) = f g^e^'^-^-^'^^e-^+^^t^c.^)) (Ag) 

3-fi 

The brackets in the average appearing in the right-hand side denote the sampling with respect to the thermalized 
initial state. We know that in this initial state 

[ 1 ^dgsinNsin^/f) if£,j>l 
(c](0) Cj (0)) = I f f^sin^-l^sin^-l)^ 1 ' \ll,j<Q (A9) 
I, otherwise 

It is convenient to rewrite 

(4(*M*)> = (4(*)Cn(*)>(l) + (4(*K(*)>(2) (A10) 

where 



a h a y 

ikm-ikn-ie^+i^t £ e+i kj (c | (Q)c . (Q)) 

4(*)c»(t)>(2) = f ^*^e <fc '' n - tt "- fe »' t + fe *' £ e- ifc '^'< C ](0) Cj -(0)> 



Let us focus, say, on (cj 7l (t)c„(t))( 2 ): 



m - ifen - iEfc ,* +fe(! t £ e - tt «-H« <c t (0) c (0) ) 

£ e- lfc ' £+4 ^ - f dg sin[^] sin[jg]/ ( 

; i 



2vr 2vr 



Now, we make use of the following identities, valid as 8 — > + , in terms of distributions: 



hence we get, the limit 6 — > + being understood, that 



— 7T 



where we have introduced the notation 

2 cos k — cos q — z<5sign(fc) 2 cos k' — cos q + £<5sign(fc') 



(All) 



(A12) 



Vrin«e ±< ^= lim ~ Slng - - (A13) 

<J— >o+ 2 cos A; — cos q =p zosign(fc) 



(cl{t)c n (t)) {2) = 2_ ^ e ^-^H^t^^j.,) (AW) 

' Z7T Z7T 



/, 2 /' 7r , „/ 9 > 1 sino 1 sino 
F(k,k') = - dqfW- ±— — — — — — A15 
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We now focus on the q integral that appears in F(k, k') : 
F(k,k') 



1 f* dq sin 2 g/f } 



2 J_ n 2ir (cos k — cos q — i<5sign(fc))(cos k' — cos q + iSsiga(k')) 
1 



cos A: — cos A:' — i(5(sign(fc) + sign(fc')) J 2ni 

1 1 



(A16) 



1 + z 2 — 2z(cos k — i<5sign(fc)) 1 + z 2 — 2z(cos k' + iSsigsi(k')) 



where we have set z — e lq and the z-integral runs counter-clockwise around the unit circle. Explicitly carrying out 
the z integral leads to 



F{k,k') = 



1 



2 cos k — cos k' — i#(sign(fc) + sign(fc')) 



ffl sin fc' + /^ sinfc 



(2) 



(A17) 



We are thus left with evaluating 



(4„(f)c„(i)) (2 ) 



w dfcdfc^ 
_ — 2tt ~2~k 



ik' m—ikn—ie k / t-\-iekt ^_ 



l 



2 cos k — cos k' — i<5(sign(fc) + sign(fc')) 



4 2) sinfe'+Z^ sinfc 



(A18) 



In order to extract the long time behavior of (|A18p . we change variables from fc to u = (k — k')t sinfc' and we expand 
in powers of 1/t keeping only the leading order at fixed m and n: 



lim (c^(t)c„(t)> (2 ) 



v dfc' Jt , 

n dfc' 



a/c „ifc'(m-n) f (2) • 
e Jfc' ' 



dw 



2n —u — isign(fc') 



2tt 



(A19) 



A similar reasoning leads to 



Hm <4(t)c n (t)} (1) = /* ^ e *(— »)/f^(fc) 



2tt 



(A20) 



Hence, (|A19[) together with (|A20[) are exactly the announced result (|A5IA6p . 



Appendix B: Green's function 
1. Green's functions for the evolution 



In order to determine Gij(t,t'), we introduce the Green's function Gij(t,t') using the Heisenberg picture, 

cj(t) = ]T G in (t, 0)c„(0), c}(t) = £ G^(t, 0)4(0) 



(Bl) 



which leads to Gij(t,t') = T, m ,n G *n(t,0)G* mj (t' ,Q)(c n (0)c< rn (0)) . In the latter formula, we take as the initial state the 
stationary-state. The Green's function G in the presence of defects can be determined by a variety of methods. We 
choose to use a path-integral formulation. Going to Grassmann fields, the Fourier transform of Gij(t,t') is given by 



Gij(w) = (ci(u)cj(w)) = 
where the action S = Sq + SS is as follows 

So[c,c] = 



JVcDc Ci (w)cj(w)e- S [ e 
j VcDce-s^ 



(-iui - h)a(uj)ci(uj) - ^CiCi+i - ^CiCj-i 



(B2) 



(B3) 
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and 

5S[c,c} = -^5hiC iCi (B4) 

and the magnetic field is Shi = 6h+6i e/% + 6h-6i^/2- We introduce gij(ui), the Green's function in the absence of 
defects as given by So- We start with the definition of the Green's function: 



(a{ui)cj{uj)) = - 



6 -Z[fj,ri] 



(B5) 



77— ?7— 



Z[0, 0] Sri, Sfk 

where r\ and f] are introduced through the partition function: 

Z[f),r)] = J V^ce- Sol5 ' c ^i.^ 5h ^ + ^ +5 ^ (B6) 

Then, we use the representation 



„y\ Shi 



J ]^P0 J P(/ )j -e~ E ^-(^^ + \/^ ( ^ Cj+Ej ^ ) ) (B7) 



giving the partition function as 

ZlViV] = J V^D^tT^ 3 ^^ 3 ^ 3 J Vtt)ce~ Sa ^^ + ^ 3 ^ 3 ^ 3+ *^^ 3 ^ 3 ^ + ( ri3+ ^ 3 ^ 



V4>V<t>e- & ^ (det g l3 )e E -.^ fe+V^Osy (% + V^7 3 ) ( B8 ) 



where we have defined 7^ = Sij — y/ ShiShjgij and Cm = \/Sh m J2j 9mjVji £n = \/Sh n ginVi- We perform the 
remaining path-integral and obtain 

Z[fj,r]} = detgijdet^je^w+Zn^inl-y-XrnU (B9) 
The Green's function can be now found from the coefficient of fjir/j in the exponential that appears in (|B9|) . so that 

Gij =gij + 22 V &hndh m gi n g m j [7 _1 ] n m 

m,n 

+ Sh+Q- - Sh-g-_)g l+ g +j + Sh-(1 - Sh+g+^g^g-j 
' (1 — $h+9++)0- — Sh-g — ) — 5h+8h-g — k9h — 
5h + 5h-(g i+ g-jg + - + g,_g +j g^ + ) 



(1 — d/i+<7++)(l — 8h~g ) — 5h + Sh-g 

where we used the shorthand notation ± for ±|. In the weak magnetic defect limit of interest here we have 

Gij =gij + Shig^g^j + Sh 2 gi + g+j 

+ 5h 2 i9i-9- 3 9— + $ h lgi+g+j9++ ( B11 ) 
+ 6h 1 Sh 2 (gt+g- 3 g+- + gi-g+ 3 g-+) + ■■■ 

It will prove more convenient to resort to the Fourier transforms of these Green functions 
G(q,q',u) =g q 5 q , q , + (sh^-M + ^e"^'-^) g q g q , 

+ f ^-(shieM-Qi +8h 2 e- i «- k '>l) (sh^"^ + ^e"^"^) g q g q ,g k + . , . 
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2. Steady-state Green's frunction 



In order to calculate the steady-state Green's frunction, we start with the expression of Gij(t,t') given right after 
(IB1[) in which the initial state is the current-carrying steady-state obtained from a step temperature profile at t = — oo. 
After Ogata [l7| . if one initially prepares the spin chain with its left hand side at inverse temperature (3_ and right 
hand side at inverse temperature j3+, one has that in the steady-state 



(4(0)Cm(0)) 



T dk 



ik(m-n) rp 



(B13) 



with the occupation number F k given by We now set out to determine Q explicitly, in the weak impurity strength 
regimes. In Fourier space we have that 



Q(Q,Q',t) 

In the weak impurity limit, given that 



dq 

2tt 



G(Q,q,t)G*(Q',q,t)F q 



G{q,q',t) =e is °%, q , 



Eq S q 



1 



£q — £k 



and retaining only the leading (time- independent) behavior in the large time limit, we arrive at 

Fn — Fqi 



Q(Q,Q',t) ~F Q 5 Q , Q , - [Shie^-Q^ +Sh 2 e-^ Q - Q ^ 

+ J g (sheW-M +8h 2 e-W- k ^) (8h x S k -Q^ +5h 2 e-* k -^) 
Fq - F k Fq> — F k 



£Q ~ £Q> 



£q - Efc £Q' — £k 



(B14) 



(B15) 



(B16) 



The above formula is at the basis of our calculation of the force. 
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